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Abstract. As recently found, the distribution of sunspots is non-axisymmetric and spot group formation implies the 
existence of two persistent active longitudes separated by 180° . Here we quantitatively study the non-axisymmetry 
of sunspot occurrence. In a dynamic reference frame inferred from the differential rotation law, the raw sunspot 
data show a clear clustering around the persistent active longitudes. The differential rotation describing the 
dynamic frame is quantified in terms of the equatorial angular velocity and the differential rotation rate, which 
appear to be significantly different from those for individual sunspots. This implies that the active longitudes are 
not linked to the depth of sunspot anchoring. In order to quantify the observed effect, we introduce a measure 
of the non-axisymmetry of the sunspot distribution. The non-axisymmetric component is found to be highly 
significant, and the ratio of its strength to that of the axisymmetric one is roughly 1:10. This provides additional 
constraints for solar dynamo models. 
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1. Introduction 

The question whether sunspots appear randomly in lon- 
gitudes has been a long-standing issue since the early 
20th century. Although the existence of preferred longi- 
tudes of sunspot formation (active longitudes) has been 
suggested long ago, the question of their persistency 
was still a subject of ongoing debates (e.g., Chidambara 
1932; |Lopez Arroyo| 1961; IBalthasar~fc ScMssler 1983; 
IVitinsky et al.| 1986; IMorcT vmov fe Kitchatinovl 2004). A 
novel analysis of sunspot group data for the past 120 
years revealed the existence of two persistent active lon- 
gitudes separated by 180° < |Berdyu gina fc Usoskinl 2003, 
BU03 henceforth). In BU03 we have shown, using differ- 
ent filtering techniques, that the active longitudes are per- 
sistent on a century time scale. An important conclusion 
of our previous work is that the active longitudes are not 
fixed in any reference frame (e.g., in the Carrington sys- 
tem), but continuously migrate in longitude with a vari- 
able rate. Their migration is defined by changes of the 
mean latitude of the sunspot formation and the differential 
rotation. Neglecting this migration results in smearing of 
the active longitude pattern on time scales of more than 
one solar cycle. In contrast with the findings of BU03, 
most previous researchers assumed rotation of the active 
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longitudes with a constant rate, which explains the diver- 
sity and contradictions of the previously published results. 

The solar active longitudes and their behaviour are 
very similar to stellar activity phenomena, including 
the flip-flop cycle detected in binaries and solar-type 
stars (IBerdyugina fc Tuominen| 1998; [Berdyugina 2004; 
|Berdyugina fc Jarvinen| 2005). On the Sun the flip-flop 
phenomenon is observed as the alternation of the major 
spot activity between the opposite longitudes with a 3.7 
year cycle (BU03). The similarity between the sunspot 
distribution and the activity patterns on cool active stars 
implies that non-axisymmetry is a fundamental feature of 
the solar and stellar dynamo mechanisms. 

The persistent migrating active longitudes imply the 
existence of a non-axisymmetric component in the solar 
dynamo and provide new observational constraints for cur- 
rent solar dynamo models. Therefore, it is important to 
quantify this effect. As found by BU03, the migration of 
the active longitudes is defined by the differential rotation 
and mean latitude of sunspot formation. In the present pa- 
per we fit this model to raw sunspot data and determine 
the differential rotation of the active longitudes. Based on 
that, we introduce a dynamic reference frame and investi- 
gate the distribution of the sunspot area in this new frame. 
This allows us to account for the migration of the active 
longitudes. We reveal a double-peaked longitude distribu- 
tion of the spot area for all sunspots without any filtering 
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and also for the single, strongest spot group observed in 
each Carrington rotation. More impressively, such a dis- 
tribution is also found for the single, strongest spot group 
observed in each Carrington rotation. Finally, we intro- 
duce a measure of the non-axisymmetry of the sunspot 
distribution and estimate a relative strength of the ax- 
isymmetric and non-axisymmetric components. The dif- 
ferential rotation law obtained and the measure of the non- 
axisymmetry can be used to constrain the corresponding 
dynamo models. Our new analysis confirms the previous 
conclusions by BU03 on a new basis and dispels the doubts 
expressed by IPelt et aD (2005) that the active longitude 
separation is an artefact of the data processing. 

In this paper we analyse sunspot group data for the 
past 120 years. We use daily data on sunspot group 
locations and areas collected at the Royal Greenwich 
Observatory, the US Air Force and the National Oceanic 
and Atmospheric Administration for the years 1878-1996, 
covering 11 full solar cycles. Here we are primarily inter- 
ested in the sunspot appearance rather than in their evolu- 
tion. Accordingly, each spot was included into the analysis 
only once when it was first mentioned in the database (ei- 
ther on the day of its birth or when it appeared at the East 
limb) , all later records of the spot were ignored. Therefore, 
in contrast to earlier studies, we analyze only the sunspot 
emergence. Because of the known asymmetry between the 
Northern and Southern hemispheres, we investigate them 
separately. About 40,000 sunspots occurred during about 
1600 Carrington rotations have been considered in each 
hemisphere. 

2. A dynamic reference frame 

As suggested in BU03, the migration of the active longi- 
tudes on the Sun is defined by the differential rotation and 
by changes of the mean spot latitude. A standard model 
of the differential rotation on the Sun relates an angular 
velocity with a helio-latitude ip as follows: 

Q = n - B sin 2 r/> - C sin 4 tp, (1) 

where f2o is the equatorial angular velocity, while B and 
C describe the differential rotation rate. Here we do not 
use the higher order term C sin 4 tp, because it influences 
mostly higher latitudes where sunspots are not observed. 
Throughout the paper, we use the time step equal to the 
Carrington period, and Eq. (JJJ takes the form 

SU = n -B sin 2 (2) 

where index i denotes the i-th Carrington rotation, and 
(ip)i is the area weighted average latitude of sunspots dur- 
ing this Carrington rotation. Based on that, we introduce 
a new reference frame which describes the longitudinal 
migration of active regions with respect to the Carrington 
frame due to differential rotation: 

i 

A 4 = A + T C Yl ( fi c-%), (3) 

j=N 
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Fig. 1. The longitudinal distribution of the sunspot area 
during cycle No. 19. The vertical axis denotes the lon- 
gitude and horizontal axis the time. The upper panel 
shows the observed Carrington longitudes and the ex- 
pected migration path of the two active longitudes (shown 
by squares and triangles) given by Eq. © with B = 3.40 
and f2rj = 14.33 degday -1 . The lower panel shows the 
same plot but after the longitude correction, i.e. subtrac- 
tion of the migration path. For better visualization, each 
plot was filtered using the Lee filter (Lee 1986). 

where Aj is the expected active longitude in the i-th 
Carrington rotation, Ao is the location of the active lon- 
gitude in the A -th Carrington rotation, Tq — 25.38 days 
is the sidereal Carrington period and flc — 360°/Tc. In 
our model Q and B are the parameters (Eq. [5J that need 
to be determined from observations, while Ao is just a 
constant defining a shift in longitude of the new reference 
frame with respect to the Carrington system in the A^o-th 
Carrington rotation. 

The model predicts that for the time Tq the 
Carrington system makes one full rotation, while the 
new reference frame does Clj/Qc rotations. The difference 
in longitude accumulated between the two systems over 
many Carrington periods results in migration of the ac- 
tive longitudes, including a number of full rotations (see 
BU03). In the following, we discard the full rotations and 
consider only the circular longitude. In the new reference 
frame, the longitude of a k-th spot in the i-th Carrington 
rotation, Afcj, is found as follows: 

Ah = Xm - Aj — n360°, (4) 

where Xki is the longitude in the Carrington frame, and 
n is defined to keep Xki within the range [0°;360°]. The 
transformation to the new reference system is illustrated 
in Fig. ^ where the upper panel shows the observed 
longitude-vs-time pattern of sunspot occurrence during 
the solar cycle 19. The gray-scale of the plot represents 
the normalized area of spots occupying a given longitude. 
We note that the cycloid-like shape of the active longi- 
tude is a typical signature of the differential rotation. The 
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Fig. 2. Longitudinal distributions of sunspot occurrence 
in the Northern hemisphere for the period 1878-1996. 
a) Actual sunspots (area weighted) in the Carrington 
frame. The distribution is nearly isotropic, b) The same 
as panel a) but in the dynamic reference frame, the non- 
axisymmetry measure is T — 0.11. c) Only the position 
of one dominant spot for each Carrington rotation is con- 
sidered, r = 0.19. d) The same as in panel c) but for 
6-months averages, T = 0.43. The solid line depicts the 
best fit double Gaussian. 

normalized area of a fc-th sunspot in the i-th Carrington 
rotation is defined as 



Aki — Ski / ^ ' Sji, 



(5) 



where S is the observed area of the spot corrected for 
the projection effect, and the sum is taken over all spots 
in the given Carrington rotation. We consider the nor- 
malized area in order to suppress the dependence of the 
spot area on the phase of the solar cycle. Otherwise, the 
analysis would be dominated by large spot groups around 
cycle maxima, with the loss of information about active 
longitudes during the times of minima. The expected mi- 
gration path of the active longitude, calculated according 
to our model is shown in the upper panel of Fig. ^ by 
squares. Triangles denote the same active longitude but 
shifted by ±180°. One can see that the expected migra- 
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Fig. 3. The same as Fig. |21 but for the Southern hemi- 
sphere. The value of the non-axisymmetry T is 0.09, 0.15 
and 0.39 for the panels b), c) and d), respectively. 



tion path follows the observed concentration of sunspots. 
Accordingly, if we correct the observed Carrington lon- 
gitudes of sunspots for the expected migration, the ac- 
tive longitudes should be constant in the new coordinate 
frame. This is shown in the lower panel of Fig. The 
sunspot area is clustered around the two active longitudes 
at about 90° and 270°. One active longitude is generally 
more active than the other at a given time, and the dom- 
inant activity switches between the two active longitudes, 
indicating a flip-flop phenomenon (cf. BU03). 

The separation of the active longitudes in the dynamic 
reference frame is clearly seen in Fig. |3 which shows his- 
tograms of the area weighted sunspot occurrence in the 
corrected longitude for the whole data set in the Northern 
hemisphere. The sunspot distribution in the Carrington 
frame shows no preferred longitudes (Fig. . The same 
distribution, but in the dynamic system, shows a clear 
preference to cluster at two corrected longitudes (Fig. |2>) . 
Note that this histogram shows the longitude distribution 
of actual spot group areas in the dynamic frame (Eq. 0} , 
without any filtering, smoothing or other processing of the 
raw data. This is consistent with our earlier result (BU03) 
that the signature of the migrating active longitudes is 
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totally smeared out in the Carrington system within 1- 
2 solar cycles, while a careful account for the migration 
of the active longitudes allows us to reveal their persis- 
tence. A more pronounced double-peaked distribution of 
the sunspot area is obtained when considering not all spots 
but only the largest spot within each Carrington rotation 
(Fig. |2t). Because of the flip-flop effect the two active lon- 
gitudes are clearly revealed. We emphasize that further 
processing and averaging of the data significantly increase 
the revealed non-axisymmetry because of the suppression 
of the axisymmetric part. For instance, semiannual averag- 
ing produces a very significant double-peaked distribution 
(Fig. Hi and Fig. 5 in BU03). Very similar histo grams are 
obtained for the Southern hemisphere (Fig. The fil- 
tered distributions (panels c and d) are shown here only 
for the purpose of visualisation, and further we will deal 
only with the raw data distribution shown in Figs. [2b and 

3. Differential rotation parameters 

Let us now estimate the parameters of the differential ro- 
tation (Eq. FJJ which fit the raw sunspot group data. For 
this purpose we employ the least mean squares method. 
First, we define the deviation between the model and the 
data as 

cm = min (\ ki ; X ki - 1 80°) , (6) 

where indices i and k denote the Carrington rotation num- 
ber and the spot group within this rotation, respectively. 
The total discrepancy is then defined as 

£ = ^EE^4, (7) 

i k 

where N=J2i J2 k A k i is the number of Carrington rota- 
tions with observed sunspots. Varying the values of B and 
£7o iu Eq. @ we search for a pair which minimizes the 
discrepancy £. Note that we also vary the value of Ao to 
obtain the minimum discrepancy. 

In order to evaluate the confidence intervals for the 
best-fitting model parameters, one needs to estimate the 
uncertainties of the observed data with respect to the fit- 
ting model. From the best-fitting parameters correspond- 
ing to the £ m in we can build the distribution of the raw 
data around the expected model (Fig. FJt) . This distribu- 
tion is nearly double Gaussian with the standard deviation 
of about a = 70°. We adopt this value as a measure of the 
random scattering of the data with respect to the model. 
Then we can compute an analogue of x 2 for our model as 
follows, 

X 2 = £lo\ (8) 

Note that this is only a scaling of £, i.e. the best-fitting 
parameters remain the same. The value of x 2 is about 0.5 
per formal degree of freedom implying that, although the 
model fits the data quite well, the data are slightly auto- 
correlated (i.e. the effective number of degrees of freedom 
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Fig. 4. The distribution of the x 2 values in the space of 
the model parameters B and sidereal Q,q (Eq. [2J for the 
Northern and Southern hemispheres. The black area de- 
picts the 68% confidence region (Ax 2 = 2.3). 

is only half the number of points). The x 2 distribution in 
the space of the model parameters B and sidereal f^o is 
shown in Fig. QJ 

The best-fitting parameters are B — 3.40 ± 0.03 deg 
day" 1 and fi = 14.33 ± 0.01 deg day" 1 for the Northern 
hemisphere, and B = 3.39 ± 0.02 deg day™ 1 and £lo = 
14.31±0.01 deg day™ 1 for the Southern hemisphere, where 
the quoted errors correspond to the 90% confidence in- 
terval for one parameter of interest with Ax 2 — x 2 ^ 
min(x 2 ) = 2.71. The results for the North and South be- 
ing close to each other are nevertheless somewhat differ- 
ent. The difference is formally significant but, taking the 
uncertainties of the parameter values obtained above as a 
rough estimate, we can only mention its indicative nature. 
The distributions in Figs. and in the dynamic refer- 
ence frame were built using the above best-fitting values 
of the model parameters. The average value of B is close 
to that obtained by BU03. We note that the result from 
BU03 corresponds to a local minimum near Slo ~14.2 and 
B = 3.46 deg day -1 in Fig.fU 

When repeating the same procedure for individual 
cycles, we obtained the best-fitting parameters B vary- 
ing from 1.5 to 4 deg day -1 , and Slo from 13.7 to 14.7 
deg day -1 (see Fig. |3J). There is a general tendency that 
smaller B are paired with larger Qq. Despite the large 
spread of the best-fitting parameters for individual cycles, 
the necessity for the differential rotation is apparent, be- 
cause no good fit can be found for B — 0. Formal averaging 
over individual cycles yields the values of B = 3.1 ± 0.6 
deg day -1 and f^o = 14.2 ± 0.24 deg day -1 . The analysis 
for individual cycles could not possibly be used to check 



Usoskin, Berdyugina & Poutanen: Non-axisymmetry in sunspot distribution 



5 



T3 J 
CD 

T3 

CD o 




5-cycle 

3-cycle 

ind. cycles 



13.5 



14 



14.5 



15 



Q o (deg day" 1 ) 



Fig. 5. The area of parameters of the differential rota- 
tion in the Northern hemisphere defined using the entire 
studied interval of 1878-1996 (large cross) as well as from 
sub-intervals: individual cycles (dotted), 3-cycle intervals 
(dashed) and 5-cycle intervals (solid). 



the persistence of the active longitudes, because each cy- 
cle is fitted independently. However, extending the studied 
interval not only systematically tightens the allowed area 
in the parameter space towards the small area determined 
for the entire data set, as illustrated in Fig. EJ but also 
proves the persistency of the phenomenon. 

In Fig. we compare the present results for the dif- 
ferential rotation of the active longitudes with measure- 
ments obtained from sunspots (Balthas ar et al.l 19861 and 
surface plasma Doppler shifts (SOHQ/MDI, ISchou et al.l 
1998). Note that the SOHO/MDI data were originally ap- 
proximated by Eq. QJ. As seen from Fig. EI at all latitudes 
the active longitudes rotate significantly slower than indi- 
vidual sunspots, which is in agreement with the previous 
finding by BU03. This implies that the active longitudes 
are not linked to the depth where developed sunspots are 
anchored. The difference in the rotation of the active lon- 
gitudes and individual spots is about 0.2degday _1 which 
produces a lag of about 2 full rotations during a solar cycle 
as reported by BU03. A comparison with the SOHO/MDI 
measurements is more difficult because of the different 
shapes of the used approximations. It appears that the 
active longitudes rotate faster than the local plasma at 
low latitudes (tp < 18°) and vice-versa at higher latitudes. 



4. A measure of the non-axisymmetry 

The existence of the active longitudes implies a long-term 
asymmetry in the sunspot longitudinal distribution which 
is related to a non-axisymmetric component of the dy- 
namo mechanism. In order to quantify this we introduce 
a measure of the non-axisymmetry in the following way. 
Let us choose the value of Ao so that the two active lon- 
gitudes correspond to the corrected longitudes of 0° and 
180° (see Fig. Eh). Depending on the value of the corrected 
longitude A^, the sunspot with the normalized area Aki 
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Fig. 6. The sidereal differential rotation of the ac- 
tive longitudes determined in this work compared with 
that obtained using surface Doppler shifts (SOHO/MDI, 
ISchou et al.l 19981 and sunspots (Baltha sar et al.l 1986). 

contributes to either of the two numbers as follows. 
A x = J2 Aki ' if 1**1 < 45 ° or l*« - 180 °l < 45 °' 



A 2 = /J , otherwise, 



(9) 



where the summation is taken over all spots in all 
Carrington rotations. Thus, N\ and A 2 represent the num- 
ber of spots, weighted by their area, which are close to 
and far from the active longitudes, respectively. The non- 
axisymmetry is then defined as 

Ai - A 2 

r = — -. (10) 

Ar + A 2 v ' 

It can take values from to 1, where corresponds to 
the axisymmetric sunspot distribution and 1 to the non- 
axisymmetric distribution. 

The non-axisymmetry of the area-weighted sunspot 
distribution for the entire studied series shown in 
Figs. Eb and is 0.11 and 0.09 for the Northern and 
Southern hemispheres, respectively. Without normalizing 
the spot area (Eq. [5}, the distributions indicate the non- 
axisymmetry of 0.07. We conclude, therefore, that the 
strength of the nonaxisymmetric component is roughly 
1:10 of that of the axisymmetric one as observed in the 
sunspot distribution for the last 120 years. For individual 
cycles r takes values from 0.07 to 0.3. We note that at 
the best-fitting parameters minimizing discrepancy £ (see 
Fig. 0}, r does not necessarily reach the maximum value, 
being, however, rather close to it. The dependence of T on 
the parameter B is shown in Fig. (parameter flo is now 
chosen to provide the maximum possible T). One can see 
that the relation has a single peak at about B = 3.42-3.43 
deg day -1 (close to the best-fitting value B = 3.39-3.40 
deg day -1 minimizing £) and decreases when deviating 
from it. If the longitudes are not corrected for the migra- 
tion (i.e., B = 0), the non-axisymmetry is T =0.02-0.03 
(see Fig. |7Ji) which is consistent with the null- hypothesis 
of the axisymmetric distribution (see below). This shows 
again that neglecting the differential rotation results in 
complete smearing of the pattern. 
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Fig. 7. The measure of the non-axisymmetry T as a 
function of the differential rotation parameter B for the 
Northern (left panel) and Southern (right panel) hemi- 
spheres. The values of B = 3.42 and 3.43 deg day -1 as 
denoted by the vertical lines, correspond to the maximum 
r in the Northern and Southern hemispheres, respectively. 
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Fig. 8. The distribution of the non-axisymmetry T calcu- 
lated for 5,000 sets of random-phase sunspot occurrence. 



In order to check the significance of the obtained non- 
axisymmetry measure T we have performed the following 
test. For each Carrington rotation we have randomly per- 
muted all the sunspots, i.e., a new random Carrington 
longitude has been ascribed to each actually observed 
sunspot while keeping its area. Then the value of T has 
been calculated as described above. We have computed 
the non-axisymmetry T for 5,000 sets of such random- 
phase sunspot occurrence. The distribution of the value 
of r shown in Fig. |H1 depicts a nearly Poisson distribu- 
tion. According to this simulation, the hypothesis of ro- 
tation of the active longitudes with a fixed rate (giving 
T = 0.02 - 0.03; see Figs. Hi, Et, and|7Ji,b) cannot be dis- 
tinguished from the null hypothesis of the axisymmetric 
sunspot distribution. On the other hand, the probabil- 
ity to obtain T — 0.11 (0.09) for an axisymmetric distri- 
bution is < 10~ 6 (w 10~ 5 ) for the Northern (Southern) 
hemispheres. This means that the non-axisymmetric com- 
ponent does really exist in the raw data of sunspot oc- 
currence and at the very high significance level. For the 
case when only the dominant spot is considered (Figs. [5]: 
andH:) T = 0.15 - 0.19, implying that nearly 60% of 
the major spots appear in the vicinity of the active longi- 
tudes. Moreover, the non-axisymmetric mode dominates 
the semiannually averaged sunspot occurrence (Figs. 
andEJi). 

5. Conclusions 

Active longitudes on the Sun were commonly expected to 
rotate with a constant rate. This a priori assumption led 
previous researchers to the conclusion that the active lon- 
gitudes, if exist, are not stable in their appearance. The 
recent analysis of the sunspot distribution revealed how- 
ever that there are two persistent active longitudes which 
migrate according to the differential rotation law (BU03). 
In this paper we have confirm this finding and determined 
the parameters of the differential rotation law affecting 
the active longitude migration. We emphasize that in the 
present paper we employ a different approach compared 
to that used in BU03, where the data were analysed to 
reveal the underlying regularities. Here we fit a theoreti- 



cal model of the differential rotation to the raw data with- 
out any pre-processing or filtering of the latter. Previously 
(BU03) we investigated only the non-axisymmetric part of 
the sunspot distribution while efficiently suppressing the 
axisymmetric one. This was criticized bv lPelt et aLl f2005') 
who claimed that the active longitude pattern found in 
BU03 is an artefact of the used method. The present anal- 
ysis, which is based on the raw observed sunspot areas, an- 
swers their criticism and confirms that the phenomenon of 
the persistent active longitudes is real. 

The main conclusions obtained in the present paper 
are listed below. 

— The raw sunspot occurrence data for the last 120 years 
confirm the existence of two preferred longitudes 180° 
apart, which migrate in any fixed rotation frame but 
are persistent throughout the entire period studied in 
the dynamic frame defined by Eq. Q. 

— The differential rotation of the active longitudes is 
well approximated by Eq. (2J with B = 3.40 ± 0.03 
(3.39 ± 0.02) deg day" 1 and sidereal n = 14.33 ± 
0.01 (14.31±0.01) deg day" 1 for Northern (Southern) 
hemispheres, respectively. This is significantly differ- 
ent from the differential rotation of individual spots. 
This implies that the depth at which sunspots are 
formed (and affected by the non-axisymmetric compo- 
nent of the field) is different from that where developed 
sunspots are anchored. 

— The measure of the non-axisymmetry T (Eq. I10|) for 
the raw data is found to be 0.11 (0.09) for the Northern 
(Southern) hemispheres. These correspond to the sig- 
nificance < 10 -6 (10 -5 ), respectively. The strength of 
the non-axisymmetric dynamo component is roughly 
1/10 of that of the axisymmetric one. Appropriate 
smoothing, filtering or processing of the data sup- 
presses the axisymmetric component and increases the 
formal measure of the non-axisymmetry. 

The above findings imply the existence of a weak, 
but persistent, non-axisymetric dynamo component with 
phase migration. The relative strength of this mode and 
the law of its rotation obtained here provide new con- 
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straints for the development of solar dynamo models as 
was recently undertaken by, e.g.. lMossl Cl999. 2004, 2005). 
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